Isolation, bioassay and 3D-QSAR analysis of 8-isopentenyl flavonoids from Epimedium sagittatum maxim. as PDE5A inhibitors

Background As known, inhibition of phosphodiesterase 5 (PDE5) has the therapeutic effect on male erectile dysfunction (ED), and the processed folium of Epimedium sagittatum Maxim. (PFES) characterized by 8-isopentenyl flavonoids is a famous herb for treating ED. However, the main flavonoids inhibitory activities, structure–activity relationship (SAR) and signaling pathway have been not systematically studied so that its pharmacodynamic mechanism is unclear. Methods We aimed to initially reveal the PFES efficacy mechanism for treating ED. For the first time, 6 main 8-isopentenyl flavonoids (1–6) from PFES were isolated and identified. Then based on HPLC detection, we proposed a novel method to screen inhibitors among them. We further analyze the three-dimensional quantitative structure–activity relationship (3D-QSAR) for those inhibitors. Results The results were verified by cellular effects of the screened flavonoids. Among 6 compounds, Icariin: (1), 2-Oʹʹrhamnosylicaridide II (2) and Baohuoside I (3) were identified with significant activities (IC50 = 8.275, 3.233, 5.473 μM). Then 3D-QSAR studies showed that the replacement of C8 with bulky steric groups as isopentenyl, C3 with positive charge groups and C4' with a hydrogen bond acceptor substituent could increase inhibitory effects. In contrast, the substitution of C7 with bulky steric groups or hydrophilic groups tended to decrease the efficacies. And compounds 1, 2, 3 could increase cGMP level and decrease cytoplasmic Ca2+ of rat corpus cavernosum smooth muscle cells (CCSMCs)by activating PKG. Conclusion 8-isopentenyl flavonoids could be the main pharmacodynamic substances of PFES in the treatment for ED, and some had significant PDE5A1 inhibitory activities so as to activate cGMP/PKG/Ca2+ signaling pathway in CCSMCs, that was related to the substituents at the key sites such as C8, C3, C4ʹ and C7 in the characteristic compounds. Supplementary Information The online version contains supplementary material available at 10.1186/s13020-022-00705-5.


Introduction
Erectile dysfunction (ED), defined as the inability of the penis to achieve or maintain sufficient erections, is one of the most common male sexual disorders [1]. Phosphodiesterase 5 (PDE5) is highly expressed in corpus cavernosum and has become a noteworthy target for drug screening in ED [2,3]. The PDE5 inhibitors relax the penile blood vessel by blocking the hydrolyzation of the second messenger cGMP. Clinically, admission of PDE5 inhibitors is the main treatment for ED. However, in specific cases, drugs containing PDE5 inhibitors were not applicable to patients, otherwise disease progression might happen, such as sildenafil could not be prescribed to patients with severe cardiovascular diseases [4]. In addition, the established side effects of PDE5 inhibitors in the market include headache, dyspepsia, vomiting, diarrhea, and vision disorders [5]. Compared to synthetic compounds, naturally derived compounds, becoming an important resource for new drugs, might exhibit higher safety value, unique pharmacological properties, and less side effects.
The processed folium of Epimedium sagittatum (Sieb. et Zucc) Maxim. (PFES), Yinyanghuo in Chinese, which is soaked in mutton oil and fried for use, is one of the most commonly used traditional Chinese medicines to reinforce the kidney Yang. In modern pharmacological research, the PFES showed ability to strengthen sexual function and treat ED [6][7][8]. The underlying mechanism was proposed that it might relax the corpus cavernosum smooth muscle through activation of multitargets on NO/PDE5/cGMP signaling pathway [9]. The characteristic chemical components of PFES are8-isopentenyl flavonoids, of which 6 compounds are high in content, including icariin, baohuoside I, epimedin A, B, C and 2-O"-rhamnosylicaridide II [10]. However, until now, it was only reported that icariin exhibited an IC 50 value of 5.9 μM against PDE5A1 in vitro [11], the 6 components have not been separated from PFES and the PDE5 inhibitory potencies have not been scrutinized. Moreover, there is few SAR study on flavonoids containing bulky substituents like isopentenyl at C8position as PDE5 inhibitors.
Compared with the detection for cGMP level using isotope labeling or fluorescent labeling, HPLC is simpler and easier to operate. Because of the fact that HPLC could detect the reduction of cGMP caused by PED5 hydrolysis, the combination of enzyme reaction activity and HPLC could be applied to analyze the activity of 8-isopentenyl flavonoids from PFES. Although recent studies presented various computational methods for drug design and discovery, 3D-QSAR analysis is still one of the most prominent and well-developed approaches, especially in predicting the activities of novel compounds and the optimization of lead compounds [12]. Moreover, based on PDE5 inhibitory mechanism, the effect and mechanism study on ED could be evaluated using cavernous smooth muscle cells (CSMCs) cellular assay in vitro. In this study, for the first time, 6 main 8-isopentenyl flavonoids, Icariin: (1), 2-O"-rhamnosylicaridide II (2), Baohuoside I (3), Epimedin A, B, C (4-6) were separated from the PFES and identified by 1 H NMR, 13 C NMR and MS, and these compounds were detected PDE5A1 inhibitory activity with the novel HPLC assay. We also adapted a series of PDE5 inhibitors from different reports [11,13] to firstly establish CoMFA and CoMSIA models, aiming at finding the relationship between the biological activity and the key structural factors of 8-isopentenyl flavonoids from the models. Furthermore, The SAR results were verified by the cellular effects of the cGMP level and cytoplasmic Ca 2+ in CSMCs, and the mechanism on PKG/Ca 2+ signaling pathway was also studied.

Plant material
The leaves of refined mutton fat-PFES were purchased and identified by Dr. Jian Xu in Jing Brand Co., Ltd., Hubei Province, China, in May 2019. A specimen voucher (20190501) had been stored in the Hubei Provincial Key Laboratory for Quality and Safety of Traditional Chinese.
The HPLC system consisted of the Agilent 1200 instrument equipped with a quaternary pump (G1311A, CA, USA), an autosampler (G1329A), a Unitary C18 5 μm reversed phase analytical column (250 × 4.6 mm, Wenling, China), a thermostatted column compartment (G1316A), a diode ultraviolet detector (G1314B) operated at 254 nm, and the analytical software for data acquisition (ChemStation B.04.01, CA, USA). NMR spectra were recorded on Bruker 600 NMR spectrometer, and TMS was used as internal standard. ESI-MS data were obtained on Agilent Series 1100 SL mass spectrometer. NS4000 System with C18 20 × 250 mm 10 μm preparative columns. Column chromatography (CC) was performed with silica gel (200−300 mesh, Qingdao Marine Chemical Inc., Qingdao, China). Sephadex LH-20 (Pharmacia Biotech, USA) was used for column chromatography. All solvents used for HPLC were HPLC grade solvents and all solvents used for extraction and chromatography separation were analytical grade solvents.

Chromatographic conditions
50 μL analytes were injected and separated on the C18 5 μm column with a flow rate of 1 mL/min at 30 °C column temperature. The mobile phase consisted of methanol and 0.05 mol/L potassium dihydrogen phosphate (10:90, V/V). The absorption was detected at 254 nm wavelength.

PDE5A1 inhibition assays
All components for assay were diluted in assay buffer which contained 100 mM MgCl 2 and 50 mM Tris-HCl (pH 8.0). The compounds were dissolved in DMSO and then further diluted to the desired concentration with assay buffer. 30 μL PDE5A1 solution and 30 μL compound solution were mixed at room temperature and incubated for 5 min, followed by addition of with 30 μL cGMP solution, and the reaction was carried out at 35 °C for 90 min. After the assay, each vial was incubated in 100 °C boiling water for 5 min to inactivate the PDE5A1 activity, and then was cooled to be detected with HPLC. No PDE5A1 was used as the positive control, and DMSO was used as the negative control. Percentage inhibition was calculated using the peak area of HPLC with the following formula.
For dose-response experiments, the PDE5A1 inhibitors were diluted to 7 concentrations series with fourfold gradient. Two individual operations were performed on different dates as the repetition.

Precision, stability and reproducibility assessment
The precision of HPLC was determined by spiking blank samples with 7.5 mg/mL cGMP (n = 6). To demonstrate the analytes stability, the analytes were placed at room temperature for 0 ~ 24 h after PDE5A1 inhibition by 10 nM sildenafil, then were injected into the analytical column for analysis. The PDE5A1 inhibition analysis by 10 nM sildenafil was repeated six times to verify the reproducibility of the assay analytical performance. The relative standard deviation (RSD) was calculated to access the analytical performance.

SAR dataset
A dataset comprising 17 PDE5A1 inhibitors was used as model dataset. The potencies of compounds 1-6 and 16 were obtained by our in vitro activity assay combined with HPLC. The remains were adapted from two literature [11,13]. Although these compounds were taken from different sources, the compound 16 [11,13] showed similar values of inhibitory concentration (IC 50 ) with ours, 0.075, 0.010 and 0.014 (our result) μM, respectively. The structures of the molecules and their biological data were given in Table 1. IC 50 of the molecules were converted into corresponding pIC 50 (-logIC 50 ). The total set of inhibitors was randomly divided into a training set (14 compounds) for generating 3D-QSAR model and a test set (3 compounds) for validating the quality of the model. The biological activity of test set molecules contains a range of low, moderate and high activity.

Molecular modeling and 3D-QSAR analyses
All molecular modeling and calculations were performed by the Sybyl program package [14]. The structures of all compounds were constructed using the Sketch Molecule module. Structural energy minimization was carried out using the standard Tripos molecular mechanics force field and Gasteiger-Huckle charge, with the convergence criterion set at 0.005 kcal/(Å·mol) and the max iterations for the minimization set to 10,000.

Molecular alignment
Because the quality and the predictive ability of the 3D-QSAR model rely directly on the alignment rules, we aligned the active conformations of the compounds derived by energy minimization previously. Compound 16 was chosen as the template molecule to fit the remaining training set and test set compounds by using the Database Align function in Sybyl [14]. Alignment of compounds in the training set was shown in Fig. 1.

CoMFA
The 3D cubic lattice region was created automatically and the grid spacing was set at 1.7 Å in the x, y, and z directions. A sp 3 carbon atom with + 1 net charge was taken as the probe atom and steric and electrostatic interactions were calculated by using the Tripos force field. A cutoff of 30 kcal/mol was adopted in calculating energies. With standard options for scaling of variables, the regression analysis was carried out using   partial least squares (PLS) method [15]. The crossvalidation with Leave-One-Out (LOO) [16] option was carried out to obtain the optimum number of components (ONC).The final model (non-cross-validated conventional analysis) was developed with the optimum number of components. Region focusing is an advanced strategy that improves prediction ability of the initial built CoMFA model. CoMFA region focusing method refines a model by increasing the weights of the lattice points. StDev × coefficient values were used as weights.

CoMSIA
Although CoMSIA is similar to CoMFA, CoMSIA not only computes the steric and electrostatic fields, but also calculates additional hydrophobic, hydrogen-bond donor, and hydrogen-bond acceptor fields. CoMSIA descriptors were calculated with the same lattice box, grid spacing and probe atom as that used for the CoMFA calculations. In CoMSIA analysis, the five different descriptor fields are not totally independent from each other. Given that the size of dataset, we decided to use three of five fields as the descriptors.

Isolation and culture of rat corpus cavernosum smooth muscle cells (CCSMCs)
Male Sprague-Dawley (SD) rats (180-220 g) were purchased from Shanghai SLAC Laboratory Animal Co., Ltd. This experiment protocol was approved by the standardizing laboratory animal ethical review of Shanghai University of Traditional Chinese Medicine. After removal of the corpus spongiosum, subcutaneous tissue and deep fascia of penises of SD rats, the corpus cavernosum tissue was cut into small pieces of less than 2mm 3 and digested with 0.5% collagenase I and 20% fetal bovine serum (PFES) in DMEM for 4 h. The cell suspensions were collected by mesh filtration and centrifugation, and cultured in DMEM containing 20% PFES and 1 × antibiotics at 37 °C in a humidified atmosphere with 5% CO 2 . After 24 h culture, spindle CCSMCs attached to the culture dish. The medium was replaced with DMEM containing 10% PFES when CCSMCs reached 70% confluence.

CCK-8 assay
A CCK-8 assay with 5,000 cells per well and five replicates in 96-plate was used to determine the toxic effects of 1 ~ 3 of 2.5 ~ 10 μM after 48 h. The inhibition rate was determined.

cGMP assay
CCSMCs were seeded in 6-well plate at a density of 100,000/well for 24 h, and then were treated with 10 μM 1 ~ 3 for 48 h. The cells were lysed by 0.1 M hydrochloric acid with 200 μl per well for 20 min incubation, then the supernatant was collected and determined the cGMP content by ELISA.

Fluorescent staining of calcium
CCSMCs were cultured in 6-well plate at a density of 100,000/well for 48 h. Then cells were stained by Fluo-4AM for 30 min prior to 10 mΜ 1 ~ 3 exposure for 30 min in the presence or absence of 10 μM Rp-8-Br-cGMPS (the protein kinase G inhibitor). After PBS rinse, the cytosolic calcium was observed by fluorescence microscope.  Tables 2 and 3, and the data were consistent to the reported literature [17]. It was elucidated as icariin.  Tables 2 and 3, and the data were consistent to the reported literature [18]. It was elucidated as2-O"-rhamnosylicaridide II.  Tables 2 and 3, and the data were consistent to the reported literature [19]. It was elucidated as baohuoside I.

Results and discussion
Compound 4. Yellow amorphous powder. ESI-MS: m/z = 861.21 [M + Na] + . The 1 H and 13 C NMR data were generalized in Tables 2 and 3, and the data were consistent to the reported literature [20]. It was elucidated as epimedin A.
Compound 5. Yellow amorphous powder. ESI-MS: m/z = 809.11 [M + H] + , The 1 H and 13 C NMR data were generalized in Tables 2 and 3, and the data were consistent  Tables 2 and 3, and the data were consistent to the reported literature [20]. It was elucidated as epimedin C.

Development of a screening assay for identifying PDE5A1 inhibitors
The development of PDE5 inhibitor, an important drug target for ED, has attracted much attention in the field. The common assays for screening PDE5 inhibitors use labeled cGMP as substrate, among which, 3 [H] is the main isobaric tags [21][22][23]. Other detection methods using fluorescence detection, such as fluorescence resonance energy transfer (FRET) [24] and fluorescence polarization (FP) [25], are also used occasionally. Although well developed, handling radioactive isotype in lab environment is rigorous while commercial fluorescence kits are expensive. Owing to these concerns, it is of practical significance to develop a novel detection method which could be popularized, economic, and easy to operate. To answer the urgent need, HPLC, as a common, simple, and effective detection technique, could be used for analyte separation and content analysis. This developed HPLC assay would be conducive for most of ordinary laboratories for PDE5 inhibitors screening, and even further be applied to studies focusing on PDE5 for other disease treatment.
PDE5A1, PDE5A2 and PDE5A3 are three isoforms of PDE5 [26,27], which have similar function and cGMP catalytic activity, and are inhibited with similar activity by sildenafil, but differ at their N-terminal regions [28]. Thus, in this present study, PDE5A1 was used as the representative PDE5 isoform in inhibitors screening from the PFES.
For developing the screening method for PDE5A1 inhibitors, the in vitro enzyme reaction and HPLC detection were combined. The decrease of cGMP peak area  was detected by HPLC according to the hydrolysis of cGMP by PDE5A1 (Fig. 2). In the range of 0.21 ~ 150 μg/ mL, cGMP standard curve fitted into well-defined linear relationship with the peak area (R 2 > 0.999). In addition, it was confirmed that the detection of cGMP was not affected by the chemical composition within buffer (Fig. 3A, B), which limited the variables in this study. The concentrations of enzyme and substrate are important parameters in the inhibitor screening assay. The results in Fig. 3C showed that the PDE5A1 activity was linear with its concentration within range of 0.025 ~ 0.30 μg/mL. PDE5A1, with concentration of 0.15 μg/ml, could potentially satisfy the requirements of acceptable signal-tonoise ratio and error at the same time (Fig. 3C). As shown in Fig. 3D, with cGMP concentration lower than 6.0 µg/ mL, linear relationship between PDE5A1 activity and its peak area was not observed due to unsaturated PDE5A1. As the PDE5A1 activity and its peak area showed high accordance with linear relationship with concentration higher than 6.0 μg/mL, 7.5 µg/mL was chosen and used as the optimized cGMP concentration for further experiments. Although the inhibition of PDE5A1 activity increased with the increase of DMSO concentration, the inhibition of PDE5A1 could be negligible when DMSO was below 0.5% (Fig. 3E). Taken together, the optimized assay conditions for PDE5A1 inhibitors screening were 0.15 µg/mL PDE5A1, 7.5 µg/mL cGMP, and 90 min reaction at 35 ℃.

Evaluation of assay quality
The quality control (QC) was evaluated for this assay. The peak area of 7.5 μg/mL cGMP was repeated 6 times. The RSD of 1.6% indicated that the instrument had good precision (Fig. 4A). 10 nM sildenafil was used to inhibit the activity of PDE5A1, then the stability of the analytes was detected by HPLC within 24 h. The stability RSD obtained was 0.46%, indicating that the analytes kept stable for 24 h at room temperature, and no obvious impurity appeared to interfere the detection of cGMP (Fig. 4B). Further, the RSD was 1.6% after six repeated PDE5A1 assay inhibited by 10 nM sildenafil, which suggested that the assay was reproducible (Fig. 4C). The signal window (SW), as the ratio of positive control signal to negative control signal, was an important parameter for evaluating the stability of screening assay. The SW mean value collected from 8 independent experiments was 4.78, which was well within the range of requirement of S/B > 3 (Fig. 4D). Lastly, sildenafil IC 50 to PDE5A1 was performed using this screening assay, and the result (IC 50 = 14.22 nM ± 1.3 nM) was similar to reported IC 50 data (5.6 nM) in other studies [28] (Fig. 4E). The QC results showed that the screening assay was stable, reproducible, reliable, and low cost for PDE5A1 inhibitors screening.

PDE5A1 activity screening of 1-6
The dose-dependent inhibition behaviors of these compounds were analyzed by PDE5A1 inhibition assay. The IC 50 of compounds 1, 2, 3, 4, 5, 6 were 8.28 μM, 3.23 μM, 5.48 μM, 408.6 μM, 357.6 μM and 682.0 μM, respectively (Fig. 5, Table 4). Among which, compounds 1, 2 and 3 showed sufficient inhibitory efficacy on PDE5A1 with compound 2 most potent (IC 50 of 3.23 μM). Our research results showed that icariin inhibited PDE5A1 with IC 50 of 8.28 μM, which was in line with other studies [11]. The high degree of consistency indicated that the novel HPLC-based detection method for PDE5A1 inhibitor was reliable and reproducible Fig. 6.  of steric and electrostatic fields were 49.3% and 50.7% in this model, respectively. For CoMSIA model, we constructed various combinations of three fields in order to gain the best model, which was shown in Table 5. The combination of steric, hydrophobic and hydrogen bond accept fields presented the optimal result. This CoMSIA model, with the optimal number of components of 7, was characterized by the following statistics: r 2 = 0.999, q 2 = 0.505, SEE = 0.042 and F value = 1663.955. The contributions of steric, hydrophobic and hydrogen bond accept fields were 21.6%, 39.3% and 39.1% in this model, respectively. The statistical parameters for the CoMFA and CoMSIA models were given in Table 6.
To validate the predictive ability of our 3D-QSAR model, biological activity of the three compounds of test set (compounds 4, 10 and 12 in Table 1) was predicted by both models. The actual activity and predicted activity data of the training set and test set were given in Table 7 for the best CoMFA and CoMSIA models. The testing results for the test set indicated that the CoMFA and CoMSIA models can be further used in new PDE5A1 inhibitors design. Simultaneously, Fig. 7 showed the correlation between the predicted pIC 50 values and the experimental ones of all the compounds for the CoMFA and CoMSIA models. The CoMFA and CoMSIA models showed comparatively satisfactory predictive capability. It was worth noting that for two models, the predicted activity values of compound 4 showed a relatively big deviation (-0.958 and -0.784), compared with its experimental value. We speculate that the lack of the data with low activity value led to the result that the QSAR model failed to identify some key features of low activity compounds.

CoMFA and CoMSIA contour maps
Contour maps, as an efficient way to explore how to increase the biological activity of the compounds studied, can provide instructive suggestions for the structural optimization of compounds. These contour maps provided us some general insight into the nature of the receptor-ligand binding region [29]. Contour maps for CoMFA and CoMSIA models were shown in Fig. 8. The most active compound 2 in the six compounds we studied, as the reference structure, was displayed in the map in aid of visualization. CoMFA The steric field contour map of CoMFA was shown in Fig. 8A. In the steric contour map of CoMFA, the larger steric features for groups near the green region or the less the ones near the yellow region were, the higher the bioactivity was. We can see that there is a large yellow polyhedron-like region around C7. This suggested that the replacement of substituent adjacent to C7 of ring A with bulky steric groups would be unfavorable for activity. This might explain why most of the compounds with a relatively minor substituent at this site, except compounds 5, 6 and 15 with a relatively bulkier group (e.g., glucose), exhibited improved potency. On the other hand, two large green polyhedrons in Fig. 8A indicate bulky substituents (e.g., isopentenyl group bound to C8 of ring A) are favorable for activity in these regions. The two substituents of ring B of compound 16 were located near the two green contours, which could partly explain the high activity of compound 16. The electrostatic field contour map of CoMFA was shown in Fig. 8B. Blue region corresponds to areas where electropositive groups were favorable, while red area indicated that electronegative groups in the position increased the activity. According to the contour maps, there are two small blue polyhedrons around the C3 of ring C, which suggested the substituent with positive charge here would increase activity. A favorable negative red contour region was found in surroundings of oxygen atom of C3 substituent of ring C in compound 2. As the experimental data exhibited, the activity of compounds 2 and 3 was higher than compound 1 because the electronegative atom in the former substituent was closer to the red polyhedron.
CoMSIA The CoMSIA steric contour map (Fig. 8C) was nearly in accordance with the field distributions of CoMFA map. A minor difference was that the yellow contour near ring B in CoMSIA steric contour map was absent. Figure 8D displayed the contour map of hydrophobic field. For the hydrophobic field, yellow and white contours highlighted areas where hydrophobic and hydrophilic properties were favored. Two white contours were showed improved inhibitory activity because of the hydrophilic group -OH of the C3 substituent present in the site. However, we could also see a yellow region around the white contour near ring A. Therefore, the length of hydrophobic or hydrophilic substituent should be selected with caution.
The CoMSIA hydrogen bond acceptor contour map was shown in Fig. 8E. In hydrogen bond acceptor contour map, the magenta and red contours identified favorable and unfavorable positions for hydrogen bond acceptors. Two large magenta polyhedrons near the C7 of ring A and the C4' of ring B were shown on the H-bond acceptor contour map, suggesting that a hydrogen bond acceptor substituent in the position may increase the activity of the compound. This was corroborated by the fact that the most active inhibitor 16 positions a substructure of sulfonyl group in the site. Similar to compound 16, compound 9, with a hydrogen bond acceptor group (-OH) at each of the two regions, displayed almost identical potency.

~ 3 increased cGMP and decreased cytosolic calcium of CCSMCs
cGMP plays a key role in regulating the relaxation of corpus cavernosum smooth muscle. In response to sexual stimulation, cGMP synthesis increases. 1 ~ 3 showed no significant cytotoxicity to CCSMCs below 10 μM (Fig. 9A), while increased the cGMP level (Fig. 9B). PKG exerts the strongest effect in cGMP-down-streamedtargeted kinase. cGMP, as a second messenger, activates PKG, which reduces cytoplasmic Ca 2+ through the cascade reactions, thereby relaxing the corpus cavernosum artery and smooth muscle and causing hemoperfusion to induce penile erection. Blockade or inhibition of PDE5 activity causes an increase in cGMP levels in the intima and smooth muscle cells of the corpus cavernosum artery, leading to vasodilation and erection. 1 ~ 3 decreased the fluorescence intensity of cytoplasmic calcium in CCSMCs, and Rp-8-Br-cGMPS blocked the down-regulated calcium levels by 1 ~ 3, suggesting 1 ~ 3 regulated cytoplasmic calcium via PKG (Fig. 9C). The results were in line with other studies, the previous studies have shown that PDE5 inhibitors promote penile erection by activating PKG [30,31]. It suggests 1 ~ 3 takes effect to CCSMCs by partly increasing cGMP, activating PKG and decreasing Ca 2+ , which may play an inhibitory role in the contraction of the corpus cavernosum smooth muscle.

Conclusion
This study was the first to isolate and identify 6 main characteristic 8-isopentenyl flavonoids from PFES and to systematically detect their inhibitory activity on PDE5. The PDE5A1 inhibition assay was carried out by the enzyme reaction and HPLC detection, and this combination was not reported in previous literature. The HPLC detection was reproducible, reasonable, simple, and effective detection tool, and was supplementary to the previous expensive detection (See Additional file 1).
The alignment-based 3D-QSAR analyses were performed to study 8-isopentenyl flavonoids structural characteristic components as potent PDE5A1 inhibitors. Both CoMFA and CoMSIA models showed satisfactory results in terms of q 2 and r 2 values. The 10 CoMSIA models with different field combination were developed. Among them, the best CoMSIA model was generated from steric, hydrophobic and hydrogen-bond acceptor fields. Besides, contour maps of the two models provided enough information about structural requirement for biological activity. Taking compound 2 as the reference structure, the replacement of C8 of the ring A with bulky steric groups, C3 of the ring C with positive charge groups and C4' of ring B with a hydrogen bond acceptor substituent could increase activities. In contrast, the substitution of C7 of ring A with bulky steric groups or hydrophilic groups tended to decrease biological activity. The proposed models could be useful for us to understand the effects of different functional groups on inhibition activities, and in future, could provide us further guidelines for designing novel PDE5A1 inhibitors with higher potency. Furthermore, it was revealed that 1-3 inhibited the contraction of the corpus cavernosum smooth muscle by inhibiting PDE5 activity, increasing cGMP level, activating PKG and decreasing Ca 2+ . This study developed a new detection method for PDE5 inhibitors. These results establish a novel HPLC assay for PDE5 inhibitors screening, provide a solid design for potent 8-isopentenyl flavonoids against PDE5, and reveal the potential mechanism by which the PFES rich in those characteristic PDE5 inhibitors is good at strengthening sexual function.